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Almost frictionless skating on ice relies on a thin layer of melted water insulating mechanically the 
blade of the skate from ice. Using the basic equations of fluid mechanics and Stefan law, we derive 
a set of two coupled equations for the thickness of the film and the length of contact, a length scale 
which cannot be taken as its value at rest. The analytical study of these equations allows to define 
a small a-dimensional parameter depending on the longitudinal coordinate which can be neglected 
everywhere except close to the contact points at the front and the end of the blade, where a boundary 
layer solution is given. This solution provides without any calculation the order of magnitude of 
the film thickness, and its dependence with respect to external parameters like the velocity and 
mass of the skater and the radius of profile and bite angle of the blade, in good agreement with the 
numerical study. Moreover this solution also shows that a lubricating water layer of macroscopic 
thickness always exists for standard values of ice skating data, contrary to what happens in the case 
of cavitation of droplets due to thermal heating (Leidenfrost effect). 

PACS numbers: 


I. FOREWORD WRITTEN BY YVES POMEAU 

The present paper reports original work which seems to me a very appropriate topic for this volume honoring 
Martine Ben Amar. Its matter is related to my first (very successful) collaboration with her. Our first joint paper jlj 
studied the growth of needle crystals in an undercooled melt. At the time it was well understood that the classical 
Ivantsov solution for the needle crystal had to be supplemented by physical effects outside of the ones in Ivantsov 
theory to yield a solution which would be both unique and pertinent for the observations. It had long been suspected 
that surface tension, through the Gibbs-Thomson curvature-dependence of the equilibrium temperature, had to be 
taken into account. But, at the time, it was more or less an article of faith. With Martine, we did show that things 
did not work as believed by many in the field: no solution exists with isotropic solid-melt surface tension. Crystal 
anisotropy had to be taken into account jl|. This advance was followed by many others thanks to Martine and her 
collaborators so that this tricky problem can now be considered as solved. 

The connection with the present work is that, as in the 1986 Europhys. Lett paper, we have to introduce Stefan 
condition for melting to understand quantitatively how a supply of heat can explain the dynamics of the solid/liquid 
interface. Ice skating relies on the same dynamics, and as the literature shows too well, it is a quite non trivial matter 
to write the appropriate boundary conditions. It is quite amazing to come back to this ’’classical”, but somewhat 
forgotten, physics after so many years. 

I wish you Martine, many more years of active and fruitful research. 


II. INTRODUCTION 

Skating is possible because of the very small friction felt by a thin blade sliding on ice, in sharp contrast with 
much larger solid-on-solid friction observed at temperatures well below the melting temperature. This has long been 
explained Q by the existence of a thin lubricating layer of liquid water. Standard wisdom was that ice melts because 
of the local increase of pressure by the weight of the skater, a pressure leading to local melting because of the pressure 
dependence of equilibrium melting temperature. But this melting is an equilibrium phenomenon, then it is unclear 
how it can describe the continuous formation of liquid water from ice, in particular because this requires a supply 
of heat. According to another explanation (jjJ the main phenomenon is melting by the heat generated by friction 
in-between ice and the sliding skate. This heat yields the energy needed to balance the latent heat, as given by Stefan 
condition on the ice surface. A condition of mechanical equilibrium expresses also how the weight of the skater is 
supported by the pressure generated inside the film, which avoids direct ice-skate contact. 

We reconsider below this question in the light of a macroscopic approach, assuming that the liquid layer in-between 
ice and skate has macroscopic thickness. Even though this layer is often referred to as a lubricating layer, it acts 
differently of a regular lubricant: lubricants like oil are chemically different of the solids facing each other. Therefore a 
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layer of molecular thickness can change the conditions of sliding. On the contrary, in the present case, the lubricating 
layer is made of the same molecules of water as solid ice, so that a layer of molecular thickness, which is about one 
angstrom for water, will be mixed with the molecules present at the surface of ice, and cannot change the ice-skate 
friction. 

Our analysis is inspired by previous work on the Leidenfrost effect [|] where a thin layer of vapor was also analyzed 
and where the balance of vertical forces (weight of the airborne Leidenfrost drop and pressure forces) plays a central 
role. A difference between the skating problem and the Leidenfrost effect is that in Leidenfrost the evaporation from 
the droplet (instead of the melting in skating) is due to the heating by the hot bottom plate although in skating the 
heat comes from viscous friction in the bulk of the lubricating layer. Here it is the longitudinal viscous friction (along 
the direction of motion) which is responsible for heating the film of water, although this longitudinal flow does not 
generate any vertical stress in the layer. The weight of the skater is actually supported by the transverse Stefan flow 
which undergoes a drop of pressure bigger than in the longitudinal direction in the skating case (this is not the case 
treated in [3], which considers a 2-D situation). 

We answer a very natural question: when a sharp edge slides on ice, how deep is the furrow due to melting arising 
itself from the heat generated by viscous friction? This introduces a length scale (the depth of the furrow) which 
turns out to be much bigger than the thickness of the liquid layer due to melting, a layer in-between ice and skate. 
To estimate this contribution of melting one needs to find out how much water is made during the time a skate passes 
over a point on the ice surface, and so how much ice has melt to make the furrow. We show that this ’’lubricating” 
layer has macroscopic thickness in concrete situations, but that it is far thinner than the depth of the furrow, an 
important remark. Furthermore, as written above, the depth of the furrow is far smaller than its length which is 
far smaller than than the radius of profile of the blade. Hence one has to solve a problem with four different length 
scales, each one with a different order of magnitude than the the three others, see equation 0 referring to the lengths 
reported in Fig|T| 

Here we present a theory using in a consistent way Stefan description of the melt/solid interface. This is not a trivial 
endeavor, because one has to write the equations for fluid mechanics in the frame of reference where the boundaries 
of this film are fixed, namely in the frame of reference of the skater. Therefore when writing Stefan condition on 
the surface of ice, we must consider the interface as a moving surface with respect to ice and the ice as moving with 
respect to the skater. It follows that the ice melting rate cannot be identified to the growth rate of the film thickness, 
as done in the recent model derived by Lozowski et al. 0-0- We show below that the missing term in their model is 
actually of prime importance. 

By combining (i) Stephan condition with (ii) the balance between the weight of the skater and the viscous pressure 
force in the liquid layer, and (iii) the relation between the volume of the trough and the melting process during the 
passage of the skater, we derive two coupled integro-differential equations relating the film thickness to the length 
of contact. All physical lengths of the problem are then deduced from the data, in particular the length of contact 
between ice and skate, a parameter which cannot be set to its value at rest as done in 0-0. The study of these 
coupled equations allows to define a small a-dimensional quantity depending on the longitudinal coordinate (in the 
direction of the motion) which can be neglected everywhere except close to the contact points at the front and the 
end of the blade. Help to this property we show that the length of contact can be deduced from the integral equation 
with a very good approximation. Then it remains to solve a single differential equation for the film thickness. From 
the analytical study of the two boundary layers at the front and at the end of the blade, we show that one can deduce 
an expression of the order of magnitude of the film thickness in these two regions in terms of the data. This provides 
without any calculation the relation between the unknown variables and the external parameters like the velocity and 
mass of the skater and the radius of profile and bite angle of the blade, in good agreement with the numerical study. 

In section m we derive the model equations for standard ice-skating conditions, namely for V-shaped blade trans¬ 
verse profile, and study in details the solution analytically and numerically with applications to hockey and inclined 
speed- skating blades. Because of the existence of a small parameter we define a critical mass below which no skating 
is possible. This critical mass is so small that one may conclude that a macroscopical film of water is always formed 
in the skating case, a conclusion not valid for cavitation of water droplet in the Leidenfrost effect Q. The case of a 
vertical speed-skating or rectangular transverse blade profile is treated in section IIVI which gives results very similar 
to the ones for V-shaped blades, although with mote complex analytical expressions. Last section is devoted to 
conclusions. 


III. EQUATIONS FOR V-SHAPED BLADES 

To get close to ice-skating conditions, we consider first a skate having a smooth surface ending toward the ice with 
a sharp edge, its longer dimension being in the direction y of the imposed speed, see figQ] Such V-shaped blades are 
the ones sharpened for hockey skaters and also figure skatings. The bottom of each blade present actually a hollow 
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FIG. 1: Schema of the blade digging into the ice surface. In (a) we set y = 0 at the front contact point and y = l at the 
back contact point where the letters t and a are placed, (b) Schema of the cross-section at a given value of y. The geometry 
is deformed in order to see the various lengths. 


separating two V-shaped edges but the skater spends most of his time on a single edge, therefore we focus here on the 
case of V-shape transverse profile on contact with the ice surface via a thin film of melted water. In this geometry the 
water layer has only a finite extend along the horizontal coordinate x perpendicular to the direction of motion, a width 
much smaller than the length of the skate. Lastly the vertical direction is associated to coordinate z. The curved 
profile of the blade and the cross-section are schematized in figsjl] where the length scales are not respected. The 
length of contact £ is defined as usual, as the longitudinal dimension of the wetted part of the blade which is situated 
between the first point of contact (where the thickness of the layer is zero) and the point of maximum penetration 
of the blade into the ice. Beyond this point the contact is assumed to stop because the bottom of the curved blade 
becomes above the surface of the through dug by the passage of the skater. This geometry also describes the case 
of speed skating (having rectangular cross-section) with inclined blade making an angle 7r/4 with the vertical axis, 
considered in subsection IIII Cl 

We neglect the fact that sharpening of the blade leaves a slightly rounded cylindrical wedge with a very small 
radius. This is a reasonable assumption because the water layer has a thickness of order of thirteen micrometers, 
although blades are sharpened with a radius of order of one micrometer. Therefore we consider a perfectly sharpened 
blade with V-shaped wedge. This hypothesis may be invalid in a very small domain (in front of the skate where the 
liquid layer is very shallow), that is not essential for our purpose which is to derive the geometry of the furrow whose 
main part has length scales much larger than the radius of curvature or the imperfections of sharpened wedges. 

We derive below the equations for the flow created by the friction forces along y which is partially evacuated laterally 
by a secondary flow flowing in the x direction (perpendicular to the plane of Fig.(a)) generating a pressure supporting 
the weight of the skater. Note that there is also some flow ejecting fluid behind the blade which is generated by the 
pressure gradient along y. This gradient is negligible with respect to the lateral one because the length of contact 
along y is much larger than the transverse dimension of the trough created by the passage of the blade. 


A. Derivation of the model with Stefan condition 

Below we write Stefan equation which links the velocity of the flow (or melting rate, normal to the ice surface) and 
the film thickness via the viscosity, the surface in contact and the latent heat of the ice. This equation set the balance 
between the heat generated by viscous friction in the liquid layer of melt water separating solid ice from the skate 
and the latent heat melting ice. This viscous friction is due to the fact that the layer has a small thickness a and 
that across this thickness the y-component of the velocity goes from V on the surface of the skate to zero on solid ice. 
This yields a viscous stress and so a source of heat in the liquid layer. For the moment we assume that all this heat 
is used to melt ice. 

Consider the slice y, y + dy. The film of water in this slice has a volume 6'(y)a(y)dy, where a(y) is the thickness of 
the film and S'(y) is the length of the wetted part of the blade, represented in FigG3-(b) as the length of the two sides 
of a triangle. The power dissipated by friction in the film of water of volume 5' (y)a(y)dy is equal to dP = rj^-S'(y)dy, 
where rj is the shear viscosity of water. In the cross section we define local coordinates (s,n), not represented in the 
figure, with — 4r < s < along the ice surface, and 0 < n < a perpendicular to the surface, n = aon the skate. 
We assume that the thickness a(y) of the film is independent of s in a given cross-section (an assumption consistent 
with the rest of the calculation) and that 


a(y) « e « £ « r, (1) 

where r is the radius of profile of the slightly curved blade which is of several meters or more (the curvature is 
intentionally exaggerated in the schema (b)). This condition implies a << S', because the bite angle 6 is typically 
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around This power dP is responsible for the melting of a volume dV = wo(y)5'(y)dy , where wo(y) is the melt 
rate. We assume that wo{y) is independent of the local coordinate s. The power dP dissipated by the shear stress is 
equal to the power required to melt the volume dV, then one has dP = Lp s dN where p a is the density of the ice (at 
the surface). It gives Stefan equation valid when the heat losses are neglected (for perfectly insulating blades) 

V 2 

a(y)\w 0 (y)\ = V— ■ (2) 

Lp s 

In the following we shall neglect the difference of density between ice and water to simplify the analysis, a minor 
approximation in regards with others (like the assumption of constant film thickness in the cross-section which will 
be assumed below), our aim being to get in fine an order of magnitude of the geometry of the trough and of the Him, 
and ultimately to define the conditions for skating with small friction. Note that the melt rate wq cannot be identified 
to da/dt , see the discussion below equation ©. 


1. Heat transfer in the film 

The equation © is a possible writing of Stefan condition. It expresses the balance between the power input coming 
from the viscous stress and the latent heat needed to melt ice. This is an approximation only, because it short cuts 
the heat transfer problem in the liquid layer. If the matter of the skate had zero heat conduction, this formula were 
exact, because the only place the heat due to viscous friction could go is to ice in order to melt it. Consider the 
opposite limit where the matter of the skate is a perfect heat conductor. In this case one has also to assume something 
about the temperature of the skate. To simplify the matter let us take it as equal to the melting temperature. In this 
case, the temperature in the liquid layer will be symmetrical with respect to the mid plane of the layer, so that the 
heat flux toward the surface of the skate will be the same as the one toward the ice. In this case only half of the heat 
generated by viscous friction will be used to melt ice. 

Moreover there is another possible effect, due to ice losses, because ice has to be heated to the melting temperature 
before to melt. As well known since Nansen [j], this has a strong influence on sliding if ice is at very low temperature: 
then the viscous friction has first to heat ice and then to melt it. Because of the quite large latent heat of ice, this is 
a negligible effect in the condition of ice-skating in arenas, where the temperature of ice is just a little below melting. 
Therefore in this case, one can assume that all the heating is used to melt ice, except of course that it is also used to 
heat the skate itself. In the case of a skate conducting heat much better than water, as argued above only half of the 
heat is used to melt ice. 

In the following we consider a blade made of perfect heat conductor, that amounts to divide by 2 the quantity of 
heat generated by the work of viscous forces and used for melting. With p s = 1, Stefan relation © becomes 

V 2 

O‘{y)\w 0 (y)\=r]—. (3) 


2. Poiseuille flow description 


Neglecting the inertia of the flow (assumed to be at low Reynolds number), we can use the analytical solution of a 
Poiseuille flow. At given y, the gradient of the pressure p in the film along the coordinate s obeys the equation 


dp d 2 u 

= V Hd 2 ' 


( 4 ) 


where u is the component of velocity in the film along s. Equation © has to be solved with the boundary conditions 
(fn)n-a/2 = u (0) = u ( a ) = «o = —V cos(0/2)^|p which is due to the geometry of the skate sunk into the furrow, 
e s k being the penetration depth of the blade at given y. The solution is 


/ s 1 dp 2 

= -na)+u 0 . 

2rj os 


( 5 ) 


The components u and w of the velocity in the cross-section are linked by the condition of incompressibility (w is 
the component in the direction n, normal to the surface of the blade). Assuming that the depth of the trough (of 
order S') and the thickness a of the film are much smaller than the length of contact £, we can neglect the gradient of 
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the pressure in the direction y of the skater motion as compared to the two other components of this gradient. The 
conservation of mass condition reduces to 


du dw 
ds dn 


Inserting the solution © in © and integrating along n in the film, we obtain 


where 


d 2 p 12yAw(y) 
ds 2 a 3 


( 6 ) 


(7) 


A w(y) = w(y, n = a) — w{y, n = 0) (8) 

results from the integral over n of the Stefan flow. The expression of A w(y) deserves some comments and requires 
to explain the meaning of the different velocities under consideration. 

Generally Stefan condition expresses a constraint on the fluid velocity normal to an interface where there is a phase 
transformation. There is a subtle point as this condition depends if it is written in the frame of reference of the 
interface or in the frame of the solid that is melting. To simplify the discussion let us consider the case, approximately 
true for ice melting, where there is no change of density from solid to liquid, as assumed here. Let us put oneself 
in the frame of reference where neither liquid water nor solid ice move, but where ice melts so that the surface of 
separation moves at speed w with respect to the frame of both ice and water (which is possible if one neglects their 
density difference). In this frame the fluid velocity normal to the surface of ice will be zero, because the phase change 
does not transport any matter. The consequence of that is that, in the frame of reference of the ice surface, which is 
moving at speed —w with respect to the frame of reference just considered, the fluid velocity near the surface must be 
also w because of the Galilean transform necessary to put the fluid in the moving frame of reference of the ice surface. 
Hence in this frame one must impose a Stefan condition to the normal velocity w of the fluid on the ice surface, or 
with our notation w(y , 0) = wo{y)- 

On the contrary, on the surface of the skate, there is no phase transformation and the boundary condition in the 
frame of reference of the skate is the regular condition (for a viscous fluid) that all components of the fluid speed 
vanish on the surface of the skate in the frame of reference where it is at rest. Finally, in the skate frame, equation 
© writes 


|Aw(y)| = w 0 (y), (9) 

a quantity which was defined as the rate of melting and cannot be identified to see relation (USD- 

Assuming the right-hand side of 0 independent of s, integrating twice over s and taking into account the boundary 
conditions (§f )o = 0 and p(6'/ 2) = 0 we obtain the drop of pressure along the blade in the cross section 

P(s) = ^(Aw/a 3 )((4) 2 - 4s 2 ), (10) 

where 8' sk is the length of the wetted part of the blade in the cross-section, equal to i e sk being given as a 

function of y in equation (ED. 

Along the vertical direction z the weight of the skater must be the opposite of the integral over s of the vertical 
component of the pressure. This yields the condition 



gs\n(9/2)Aw{y) 



Mg , 


( 11 ) 


where £ is the length of contact, M the mass of the skater sliding on a single wedge of the blade (and assuming 
that its weight is supported by a single skate) and g the acceleration of gravity. 
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Now we have to describe the formation of the trough, that provides a relation connecting wq and V . The total 
area of the melted water in the V-shaped cross-section at y is equal to ^e(y)S(y). This area results from the passage 
of the skater and the ensuing melting during the time St = y/V. Therefore 

Ay ) = \ z { y ) Kv ) = d y ' wo ( y ') s '{ y ')/ v , (12) 

or in a differential form, (after setting S'(y) = c ^g^ 2 ) an( i 'Hi/) = 2e(y) tan(0/2)) 

de 

wo{ y ) =-7^Vsm6/2, (13) 

where e(y) may be be written as the sum of two terms as e(y) = e s k(y ) + B in(e/ 2 ) ■ this f° rmu l a e sfc is the penetration 

depth of the blade at a given y in the frame of reference of the skate and si ^gj 2 ) ' s f° r the vertical height of the liquid 
layer, a layer which originates also from the melting of ice. 


3. Length of contact and film thickness 


Using equations (f3|)- (llll) - (ll3l) we can derive now the geometry and the components of the velocity in the film below 
the skate. Let us first calculate the length of contact £ using equation (fill) , that requires to express A w(y) and 
The blades of hockey and figure skates are curved along the propagation direction, in the plane (z, x), with a radius 
of profile r. This radius of profile being much larger than the length of contact (for hockey skating r is generally 
between 2m and 4m), the penetration depth of the blade writes 

2/(2£ - y ) 


f-skiy) = 


2 r 


(14) 


In a given cross-section the wetted length of the blade is S' sk = rcos(e'/ 2 ) • Defining \w s k\ = V sin (0/2)^^, which 
could be seen as the normal velocity of the blade with respect to an observer at rest (or the inverse), we have 

V 


and the flow velocity wq writes as 


\w s k(y)\ = V sin(0/2)- 


wo(y) = \w s k{y)\ + 

ay 


Introducing equations (fl5l) - (fTOl) in Stefan relation (JoJl, gives the differential equation for the him thickness a 

,£ — y da' 

' + ~dy. 


a ( sin(0/2)- 


= fc, 


(15) 


(16) 


(17) 


which has to be solved with the initial condition a = 0 ahead of the blade (at y = 0) because we have neglected the 
longitudinal squeezed flow. In (ED the length scale k is given by the relation 


k-'A 
k - 2L ’ 


(18) 


a very small length in the case of lubricating melted water him and realistic velocities (equal to 1.810 8 m for V = 
12m/s), therefore our study takes place in the limits 


fc«a«e<C£«r. 


(19) 


As found below, the him thickness is of order 15/z, the penetration depth of order 1mm, the length of contact of order 
10cm, and the radius of profile of order of a few meters. Equation ED is to be considered jointly with the set of 
equations Q and (fT5l) - (fT(jl) . together with the condition ED balancing the skater weight with the beneath pressure. 
This later set can be reduced to the integral relation, 


Mg = 


r\V sin(0/2) 
k 3 


2 £y - y 2 
r cos(0/2) 


sin(0/2) 


t-y 
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( 20 ) 


Our mathematical analysis finally yields the coupled set of equations ED-fEHD- This makes our main result and it 
will be studied in the rest of the section. 
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B. Solution of coupled equations 11171) and (1201) 

Equation ( 0 ( 71 ) fixes the value of l once the first equation 1171 ) is solved. In principle this makes a well defined 
problem. Since the equations are nonlinear, no analytical solution exists and a numerical solution must be looked at. 
This can be done by iterations by taking successive values of £ until the procedure converges: choose an arbitrary 
value £^ for the length of contact, solve the differential equation 03 which gives a solution , insert this solution 
in the r.h.s of equation ([HD and increment £ using the same procedure as before until condition COD is satisfied. 
However in the situation of ice skating, because the length scales are widely different, it makes sense to try to find 
the solution in the frame of the inequality (flUl) . 


1. Solution at leading order 

From inequality ID one may conjecture that equation HD has to be solved at leading order, by considering the 
derivative as small in the term ^ ^Defining 


«(y) 


da . rda . 

- /wo = V—/»„, 


( 21 ) 


we get the simple solution 


a (0) (y) 


sin( 0/2 ){t-y))' 


( 22 ) 


in the limit k <C 1 , where the superscript (0) is to mean that this is the solution at leading order. Such an 
approximation requires to define a small parameter ( 7 , equation (1251) 1 in terms of which the expansion of the solution 
is justified, at least in a certain domain. This solution has two obvious defects. It diverges at i = y which contradicts 
© and makes infinite the derivative ^ which was assumed to be much smaller than to derive equation (ED from 
equation ED .Therefore there is a boundary layer near ( = y (physically where the blade becomes exactly horizontal). 
The second defect of the leading order solution (1771) concerns the vicinity of y = 0 where there is a second boundary 
layer, because at the front part of the ice-skate contact a(y) must tend to zero, a limit value obviously not satisfied 
by the solution (1771) . 

However it turns out that, because the two boundary layers have a small extension in the y direction one may neglect 
their contribution to the integral in equation (1201) , a result valid for all skating configuration as discussed in subsection 
1IIIB 21 This integral which has to be calculated as well by canceling the derivative becomes independent of a(y) 

and can be computed. Using the result J Q f dy{£ — y) 4 (2£y — y 2 ) 3 = we obtain 


£(0) = (cw CMg )' (23) 

with 

(cos(6>/2 )) 3 (rjVfr 1 

sin( 6 »/ 2)) 5 ( 2 L ) 3 ' 1 } 

This expression indicates that the length of contact is almost independent on the mass of the skater (as the power 
M 1 / 11 ), depends on the velocity as U 2 / 11 and mainly depends on the radius of profile (as r 7 / 11 ). Note that ultimately 
the heat losses have a quite small effect on the water layer geometry, because it just amounts to reduce the length of 
contact by a factor 2 3 / 11 ~ 1 . 2 . 

To close the leading order solution of the coupled set CD and (ESD, one has to replace of £ by in (1771) . Note 
that by plugging into this expression of the numerical values for hochey skating, one finds that the inequalities 
©2) are satisfied. In terms of the data, the approximation that the derivative ^ is negligible in equation 122D , writes 
7 C 1 with 



( 25 ) 







because, by assuming ED. we have = ~ sin (g/ 2 1 )^_ y ) 2 ) which is of order pp in the major part of the domain 

and the term sin( 0 / 2 )-^^ is of order £/r (the bite angle 6 is a parameter which can be taken away from the discussion 
of order of magnitude because it is about 7t/2). If the inequality kr 2 -C £ 3 is satisfied, the leading order solution (1771) 
for a should be valid outside the boundary layers. This inequality is fulfilled for skating. This is not a consequence 
of k <C £ ■€ r, but follows from the order of magnitude typical of data, this ratio being of order 10~ 3 for figure, 
hockey and speed skating. However the condition 7 1 does not guarantee that the leading order solution (1771) is a 

good approximation for £, this point will be discussed in subsection IIII B 2 l where another condition, more drastic, is 
derived. 

To go further of the understanding of the role the y-dependent ratio n(y) we shall investigate the behavior of the 
solution close to the two boundary layers. 


2. Close to the boundary layers 


Close to y = £ and y = 0 local forms of equation m can be analytically derived by using standard asymptotic 
expansion of a solution of a differential equation with a small parameter. This requires to define specific scalings valid 
in regions much smaller than £. 

1) At the front of the blade, in the limit y — > 0, the solution of equation (1T71) with i.c. a(0) = 0 is a(y) = (2 ky) 1 / 2 , 
which has a square root singularity. In a certain domain one can neglect y with respect to t in (1171) . This leads to 
define, in this domain, the scaled variables y = and a = -^ with the length scales 4 = and an = , a scaling 

consistent with the use of the lubrication approximation for the liquid layer because £ <C r. In the intermediate 
domain which makes the transition between the solution a(y) = (2 ky) 1 / 2 and the exact one, equation ED takes the 
local form 


S (sin{ 0 / 2 ) + 


= 1 . 


(26) 


Its solution behaves asymptotically like 

1 

01 sin( 0 / 2 ) 

which matches the leading order solution ( 1221 ) far away from the front part of the skate in these scaled variables. 
Restoring the original variables, the solution a bL (y ) of the local equation (l26l) is undistinguishable from the exact one 
for 0 < y < £q, as illustrated by the insert of Fig. [2 Further the local solution tends to the constant value ao = jr- 
The range £0 << y « £ defines the matching domain. At larger values of y, the local solution remains constant 
whereas the exact solution a(y) remains close to a^°\y) until the second boundary layer. 

2 ) At the end of the contact, where y tends to £, the leading order approximation is not valid, as written above, 
because the solution (1221) is derived under the condition that the derivative ^ is small, whereas this derivative diverges 
at y = £ (as well as a). Therefore we have to consider that all terms of equation ed become of the same order of 
magnitude, that leads to introduce the following length scales 


a* = (fcV) 1 / 3 ; 4 = (kr 2 ) 1 / 3 - (27) 

which gives in particular the order of magnitude of the thickness in this domain. Notice that because tCr the 
length scale a* is much smaller than 4 as necessary for maintaining the validity of the lubrication approximation for 
the flow in the layer between the skate and ice. 

By introducing y = /-p 1 - and a = pp equation (1171) is transformed into 

a (sin(6/2)y - = 1, (28) 

which is in fact identical to the original one, up to the rescaling defined in ED. and the change of variable from 
£ — y to y. Therefore in this case the local solution is the same as the global one and we cannot derive any analytical 
solution or numerical one from initial conditions taken at y = 0 (y = £), because the forward solution a of the 
differential equation C51) is unstable (as well as the backward solution a(y) of ED). This result could appear without 
any interest, but it is not, because the relation (1271) is quite relevant, since it gives the order of magnitude of all the 
unknown lengths without solving any differential equation, in agreement with the exact numerical results, as reported 
next. 
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FIG. 2: (a) Film thickness a(y) and (b) penetration depth e(y). The leading order solution (blue-dashed line) is compared to 

the exact solution of equations m -ct with i.c. a( 0 ) = 0 , drawn in red-solid line which is indistinguishable from the first 
order one. 


C. Numerical results 

In the numerical study we take the following data (in MKSA units) = 10~ 3 kg/(ms), L = 3.3 10 5 W/m 3 , 
g = 10 m/s 2 and for the skater V = 12 m/s, M = 75 kg (except in the subsection devoted to the role of these 
parameters). Moreover we set respectively r = 3m and r = 25 m depending on wether we describe figure/hockey 
skating or speed skating (the radius of profile being between about 2m for figure skating, between 2m and 4m for 
hockey, and may be longer than 25?n for speed skating). In all cases we show that the assumption 7 < 1, made to 
guarantee that the parameter n(y) is small in the major part of contact, is fulfilled. More importantly we are going to 
show that the leading order solution £^ is a good approximation for £, that allows to shorten the numerical calculus 
of the solution of the two coupled equations by making a single iteration. 


1. Figure and Hockey skating 

The exact solution of the coupled equations requires to increment £ step by step, using the method outlined above: 
solve calculate the r.lr.s. of (0(71) until it is equal to its l.li.s. We find that the resulting length of contact is very 
close to . Therefore a short cut to get a solution in good agreement with the exact one is to use this leading order 
value for £ and solve 113 only once. Doing this the solution of equation ED with initial condition a( 0 ) = 0 is found 
to agree very well with the exact one (the error being less than one per thousand for the film thickness), so that they 
are indistinguishable in fig. [2] both merging in the red solid line. The numerical results obtained with the single step 
integration of (1171) ) are 


a < - 1 - > (£) = 14.83/r; e ( 1 ) (f) = 0.87mm £ {0) = 7.13cm; (29) 

in very good agreement with the exact values, 

a(£) = 14.83/x; e(£) = 0.86mm; £ = 7.10 cm, (30) 

obtained by the multiple step method outlined above. These lengths are ordered as required by equation m 
moreover we have wo <C V, as expected. 

For comparison the leading order solution is shown in blue dashed line. In (a) ( y) agrees with the exact one 

over the main part of the bottom blade, but differs in front close to y = 0 (hardly visible in this figure but detailed 
in figHI) and close to y = £, as stated previously (subsection IIIIB ll) . To illustrate the validity of the conjecture made 
(at the beginning of the previous subsection) that the derivative j// is small compared to we compare these two 
quantities in Fig. [3] In this figure the derivative term is enhanced by a factor ten to be visible. The dashed curve, 
drawn for 10 -^, shows that the derivative term is at least one order of magnitude lower than ^ (solid line), except 
in the vicinity of y = 0 and y = £, as expected. This result shows the importance of taking account of the relative 
velocity of the skater with respect to the ice surface, when writing the expression of the melting rate, as yet noted. 

In the two boundary layers the numerics agree with the analytical study presented just before. Close to the origin 
Fig. 0]shows the solution of the local equation (1261) together with the exact one. They merge including in the matching 

7 , 2 

domain which extends on a distance few times £q where £q = ^j/r = 0.032 mm. The insert displays the whole domain 
££<y <£ to show that the local solution is no more valid beyond the matching domain where it remains constant. 
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FIG. 3: (red solid line) and ^ versus y (blue dashed line) 


a,™™,, 



FIG. 4: Boundary layer description close to y = 0: the exact solution (red solid line) is compared with the local solution 
of equation m a b L (in restored variables, plotted in purple dashed line). They merge in the principal curve drawn in the 
domain 0 < y < 20£o. £o = 0.03 mm and ao = 0.77 y. The insert shows the solutions over the full domain to show that the 
matching domain (defined in the text) is much shorter than £. 


The second boundary layer, close to y = £ , is clearly visible in Fig|2]-(a) because the leading order solution a*- 0 - 1 is 
well separated from the exact one over a noticeable distance. The extension of this layer, of order £* = 0.5cm (the 
distance between the two vertical lines), is much larger than the extension of the boundary layer in front. We point 
out that the analytical study of the boundary layer presented above allows to derive straightforwardly the order of 
magnitude of the thickness in terms of the data. Equation m gives a* = 10/i in qualitative agreement with the 
exact result, a result obtained without solving the coupled set of equations (ED-®- 

To compare our theory with the one presented in references we note that they differ on several points. 

First the condition m necessary to describe the formation of the trough is lacking in this study, that leads them to 
estimate the length of contact from its value at rest, whereas the dynamical length of contact obtained here is five 
times larger. Secondly the authors of i-i identify the melting rate wo with the thickness growth rate % although 
we show that it is more than one order of magnitude larger almost everywhere below the blade. This result shows the 
importance of taking into account the relative motion of the ice with respect the skater when writing the expression 
of the melting rate w o in the frame of the skater. Lastly the authors of !-! add/subtract linearly the effects of 
the different contributions to the thickness growth rate, that leads to a single differential equation, although our 
treatment leads to a set of two coupled equations, one for the film thickness and the other for length of contact. 


2. Role of the parameters 

It could be interesting to know the optimal conditions of skating. To give an idea of the role of the various 
parameters on the film thickness, let us investigate the role of the parameters which could be changed, like the mass 
M and velocity V of the skater, and the sharpening of the blade (its radius of profile r and bite angle 9). Here we 
investigate the role of these external parameters, keeping unchanged the other ones. The result is summarized in fig. 

m 

Our curves show clearly that the film thickness and also the length of profile strongly depend on the velocity of the 
skater, moderately depends on the radius of profile and on the bite angle, although the geometry of the layer is almost 
insensitive to the skater mass. Let us now compare these numerical results with the analytical study of solution of 
sec. IIIIB 21 Introducing the bite angle dependence, the order of magnitude of the film thickness (at y — £) given by 
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FIG. 5: Film thickness a(£) and length of contact £ versus the following data: (a) skater velocity U (all other data unchanged), 
(b) radius of profile, (c) skater mass and (d) bite angle; in (b)-(d) t/=8m/s. 


(EZl) becomes, 


a* = 


/A) 2/3 ( rV 2 
y 2L > \sin (0/2) 


1/3 


(31) 


This relation is in qualitative agreement with all numerical results displayed in figE] It predicts a stronger dependence 
versus V than versus r and 9 , and independence with respect to M. Therefore we can claim that it contains the main 
physical result of the present study. 


D. Speed skating with an inclined blade 

The model derived above can be applied to the case of a speed skater strongly inclined with respect to the vertical 
axis, this angle being equal to 7r/4. A more general formulation can be derived in the general case of any angle value, 
but not presented here. For an inclined and symmetric rectangular blade , the geometry of the contact region is very 
close to the one considered for hockey or figure skating, but has to be calculated for larger radius of profile value and 
setting 9 = 7r/2. The radius of profile for speed skating is almost ten times larger than for figure/hocker skating, its 
typical value being 25 to, although the value of 9 is close to the typical values of hockey-skating for which bite angles 
are between 70 and 85 deg. 

Because the derivation of Stefan flow is the same for an inclined blade as for a V-shaped blade, it is enough to 
insert these two different parameters 9, r into the equations obtained in section cni We get the following results , 

a(£) = 29.13/z; e(£) = 1.35 mm; £ = 25.6cm, (32) 

This shows the role of the radius of profile: increasing it by a factor 8 increases by a factor 2 the depth of the through 
and the thickness of the film. The geometry of the trough is detailed in Fig. [6]where the exact solution of the coupled 
equations (I13D - (I2()|) are drawn in solid lines, the leading order solution is drawn in dashed line for comparison. In 
summary we show that speed-skaters can plough a trough twice larger and deeper than hockey or figure skaters when 
skating with an inclined blade at the same velocity, and that the film of water is also twice thicker. These results are 
completely different from the study of [g] where the film thickness is of order 1/t, whatever the angle of the blade. 


E. Critical mass 

The analytical study of the coupled equations ED -EHD has enlightened the role of the quantity k( y) defined in 
m■ We have shown that this quantity can be neglected to get the solution of the integral equation CU1) : this gives 
£( 0 ) 

in very good agreement with the exact length of contact, although k cannot be neglected to get the solution 
of the differential equation ED- because it becomes of order unity, as discussed in section IIII B 21 and illustrated in 
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FIG. 6: Solution for a speed skater making an angle 7r/4 with respect to the vertical (a) Film thickness a(y) , (b) penetration 
depth e(y), (c) velocity of the melting flow normal to the ice surface in the through. The exact numerical solution of equations 
go -m are drawn in red solid lines and the leading order solution in blue-dashed lines . 


Figs(2|a)-[3][11 There remains a question: What is the domain of validity of the approximation £ ^ £° ? In other- 
words is there a range of parameters such that the leading order solution £^ differs noticeably from the exact result? 
To answer this question we shall quantify the contribution of k in the integral equation (1201) by using the results of 
subsection IIII B 2l Let us consider the integral in the r.h.s. of (12U1) 


1 = 



(33) 


and set AX = X — X where X is for w o = Wq°^ . The condition -C 1 is the one we are looking for, because it 
means that £^°' ) is a good approximation of l. The integrand in AX is significant only in the region y ~ £ of extent £*. 
Assuming that k ~ 1 in this domain, and ignoring the role of the bite angle 9 we obtain 

0 

AX = 15(—) 4 K*( —) 3 , 
r r 


moreover we have X ~ nfs yr at leading order. In terms of the length scales, the condition <C 1 can be written as 

j « 0.25, (34) 

a condition fulfilled above (the ratio is equal to 0.03 in the case of FigEJ)- To compare with the inequality ( 1251 ) 
which is for the condition k <C 1 in the central domain, we may notice that 7 = (- f ) 3 , therefore the condition (1341) is 
more drastic than (USD- This is because for deriving (1341) we make one step more since we impose that the contribution 
of k in the boundary layer is negligible. In terms of the data, using the order of magnitude (Em for K* and equation 
m for £, the condition (El allows to define a set of critical parameters for the mass or velocity of the skater, and 
for the radius of profile of the blade. For example the inequality (134)) writes M M* (V, r) with 

M* ~2L(rfc 5 ) 1/3 (35) 

where M* has the dimension of a mass and the length k is proportional to the speed V. Using the data taken above 
for hockey skating with V = 12m/s , r = 3m, this critical mass M* is very small (0.2g). To get critical mass values 
of order one kg, the velocity should be increased by a factor 100 which would be inaccessible for skaters. This result 
shows that the relation ( 1271 ) covers all realistic skating situations. 

From the previous estimates, one can also derive the friction force on the skate of length l. This is the friction on a 
surface 16' in a fluid where the velocity gradient is of order V/a. This total force (not per unit length of the skate) is 

_ V V16' 

-T V — 

a 

where S' and a have been estimated above. A dimensionless measure of this friction is the ratio of F v to the weight 
Mg of the skater With the data taken above, this ratio is about 2.10 -4 . 

As a final remark, let us note that the geometry considered in this section, namely a narrow furrow made by melting- 
ice, solves the problem of turning, because by tilting at the right angle, the skater rests on an inclined furrow which 
can stand a priori both his/her weight and the centrifugal force generated by the turn. 
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IV. RECTANGULAR BLADE: VERTICAL SPEED SKATING 

We keep the same notations as before and assume that the skater moves on a vertical (not inclined) blade, with 
a rectangular cross-section of width <5 independent of y. The blade is sunk into a furrow of total depth e(y) = 
e s k(y) + ce(y), where a is the thickness of the film of melted water. To simplify the calculations we assume that the 
thickness is constant all along the cross section, namely that a only depends on y as above, whereas the value of a(s, y) 
could be larger vertically than horizontally where the film is squeezed by the weight of the skater. Then the furrow 
has a total depth e(y), total width S + 2 a(y), and the length of the ice surface at the abscissa y is S'(y) = S + 2 e(y). 

Equation (0 which describes the balance between the heat dissipated by viscous friction and absorbed to melt 
the ice surface is still valid, so are Poiseuille equations (g])-© and continuity condition ©, together with equations 
Q-H2D- The relation between the melting rate wo and the slope of the skate blade is deduced as above by writing 
the condition for the formation of the trough, A(y) = e(y)S = / 0 y dy'wo(y')S'(y')/V, that gives after derivation with 
respect to y, 


w 0 (y) = V 


de 


6+ 2e{y) dy 


The drop of pressure m is also the solution of 0 (with the boundary conditions as in previous section), 
integration of the pressure over the surface of contact, balanced with the weight of the skater writes here 


(36) 


The 


Mg = r/S 


, (y) j 

wo{y) ,, s dy, 

a 6 (y) 


(37) 


where K\{y) = S 2 + 6 (e 2 (y) + e{y)S). 

Finally, setting e = e s k + a in (1371) and (1M1) . the two coupled equations relating the thickness a(y) to the length of 
contact £ are 


which replaces and 



fc( 1 + 


2 £y - y 2 
Sr 


) 


which replaces COl) . with 


and 


= d «(TY Kllv) 

wo = s ^-y) + il r 

V Sr + (2 iy — y 2 ) + 2 ra(y) 


K\{y) = S 2 + 6 ( ( 


yi - y 1 


) 2 + S 


yl - y‘ 


(38) 


(39) 


(40) 


1. Length of contact at leading order 


Neglecting the film thickness and its derivative in the above equations (1551) - (HU1) gives the leading order solution 
' de(0) With e(°) = (2£(°) y _ y^/r we obtain 


i5+e(°) dy 


1 +- 


and 


°Y°\y) = k 


Sr — (2 t^y — y 2 ) 
(£(°) — y)S — 2 hr 


(41) 


and a length of contact solution of the integral relation 
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FIG. 7: Film thickness and penetration depth in the case of speed skate straight on ice (rectangular bottom blade) (a) Film 
thickness a(y) , (b) penetration depth e(y), where the exact numerical solution of equations (1381) - (13911 are drawn in red solid 
lines and the leading order solution in blue-dashed lines . 




2 £ { 0 ) y-y 2 ^ 

r 




Mg 


k 3 

Vrjd 


(42) 


which can be solved by iteration. The numerical result is £^ = 22.2cm, and a curve cd°)(y) which diverges at 
y = a® (y) as in the case of V-shaped blade studied above. 


2. Exact solution 

The exact solution may be obtained as above by incrementing , solving (1551) with £ = £^ to obtain a ^, report 
this function in the r.h.s. of (1551) until the right-hand side equals Mg. But we may also just solve equation (1551) with 
initial condition a(y) = 0 assuming that the length of contact is t = £^ because it gives quite the same result for 
typical data, as in the case of V-shaped blades. This procedure gives a film thickness and trough geometry of size, 


a{£) = 52/r; e{£) = 1.0mm; i = 22.2cm, 


(43) 


as illustrated in Fig. [7] 

Note that these numerical results are similar to those of the inclined blade. Without making a detailed study of 
what happens in the general case of a rectangular blade wedged dissymmetrical into ice with any angle, we can make 
the conjecture that the numerical values of the lengths will be of the same order as those reported in (1521) and (1451) . 


3. Estimation of a(£) 

Because the calculations are heavier in the case of rectangular blade than in the case of V-shaped blade, it could 
be tempting to get the order of magnitude of the film thickness directly from the original parameters. This can be 
done as above by looking at the scalings close to y = i. There, defining the length scales for i — y and a* for a, we 
assume that all quantities in (1551) have the same order of magnitude. Actually we may neglect the term — which is 
small with respect to t*/r, (a condition that has to be verified a posteriori), that gives = (fcr 2 ) 1 / 3 , and 

a* = (rP) 1 / 3 (44) 


where 


k = k(l + —). 

ro 

Taking the above data it gives £* = 3 cm and a* = 39/r which still qualitatively agrees with the numerical value of the 
film thickness close to y = £. We must notice that a* cannot be easily expressed in terms of the data, as it was in 
the case of V-shaped blade, because k depends on the parameter £ whose expression at leading order, equation (l42l) . 
is not simple, equation HD 
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V. SUMMARY AND CONCLUSION 

The problem of understanding how and why it is possible to skate on ice almost without friction has remained 
imperfectly understood for a long time. In this work we introduced and discussed in depth all the necessary ingredient 
of a coherent theory. In particular we emphasized that, for a curved wedge-like blade, the one used in all skating 
sports, either speed, hockey or figure skating, the geometrical parameters and the velocity of the flow are linked 
together by three rather complex relations. One describes the melting of ice by viscous friction in the thin layer of 
melted ice (Stefan condition), another one expresses that direct contact between the skate and ice is avoided because 
the pressure of the Poiseuille flow in the layer is sufficient to lift the weight of the skater above the ice surface, the 
third one is the kinetic condition which connect the volume of melted ice to the volume of the furrow dug by the 
blade. Putting together those three relations we have shown that all the unknown physical quantities can be deduced 
from a set of two coupled equations, one differential the other integral, for the film thickness a and the length of 

da 

contact l. The analytical study of both equations has revealed the existence of a parameter (k = which plays 
a key role in the derivation of the solution because it is small all along the film, except at the ends. This property 
shows in particular that the melting rate wq definitely hasn’t to be confused with the growth rate of the layer ^, an 
assumption found in the literature, because the two quantities differ by one order of magnitude in the correct theory. 
Besides this remark the fact that such a small parameter exists is of prime importance from a practical point of view, 
because it allows to get analytically the order of magnitude of the solution and inform about the effects of the input 
data. 

A natural extension of this complete theory is to allow the design of improved skates by optimizing their shape to 
lower friction and increase the grip on the ice surface when the skater makes a turn. 


VI. ACKNOWLEDGMENT 

One of us (YP) would like also to thank his grand son, Gaspard, for teaching him patiently the ways of hockey 
skating. 


[1] M. Ben Amar and Y. Pomeau, Theory of dendritic growth in a xeakly undercolled melt, Europhys Lett, 2 (1986), p. 
307-314. 

[2] O. Reynolds, Papers on mechanical and physical subjects, vol.2, p. 734-738, Cambridge University Press (1901). 

[3] F. P. Bowden and T. P. Hughes, The mechanism of sliding on ice and snow, Proc. R. Soc. London, Ser. A, 172 (1939), p. 
280-298. 

[4] Y. Pomeau,M. Le Berre, F. Celestini and T. Frisch, The Leidenfrost effect: From quasi-spherical droplets to puddles, C.R. 
Mecanique 340(2012) p. 867-881. 

[5] E. P. Lozowski and K. Szilder, Derivation and new analysis of a hydrodynamic model of speed skate ice friction. Int J 
Offshore Polar Eng 23 (2013) p.104-111 

[6] E. P. Lozowski, K. Szilder and S. Maw, Sports Engeneering, DOI 10.1007/sl2283-013-0141-z, published online October 5, 
2013. 

[7] F. Nansen, Farther north, 1 (1898), George Newnes, London. 

[8] A.J. Fowler and A. Bejan, Contact melting during sliding on ice, Int. J. Heat Mass Transfer; 36 (1993) p. 1171-1179. 

[9] L. D.Landau and E.M. Lifshitz, section 9 of Theory of Elasticity, Pergamon Press, Oxford (1970). B. Audoly and Y. 
Pomeau, section 14.2 of Elasticity and Geometry, Oxford University Press, Oxford (2010). 

[10] S. Mora et al. Solid drops: large capillary deformations of immersed elastic rods, Phys. Rev. Lett. Ill , 114301(2013) 



